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ABSTRACT 

We investigate the dependence of the Fanaroff- Riley (FR) I/II dichotomy of ra- 
dio galaxies on their luminosities and redshifts. Because of a very strong redshift- 
luminosity correlation (Malmquist bias) in a flux-limited sample, any redshift- 
dependent effect could appear as a luminosity related effect and vice versa. A question 
could then arise - do all the morphological differences seen in the two classes (FR I and 
II types) of sources, usually attributed to the differences in their luminosities, could 
these all as well be a result of mainly a cosmological evolutionary effect (e.g., due to 
the changing ambient density) with cosmic epoch? Even a sharp break in luminosity, 
seen among the two classes, could after all reflect a rather critical ambient density 
value. A doubt on these Hues does not seem to have been raised in past and things 
have never been examined keeping this particular aspect in mind. We want to ascer- 
tain the customary prevalent view in the literature that the systematic differences in 
the two broad morphology types of FR I and II radio galaxies are indeed due to the 
differences in their luminosities, and not a manifestation of an evolutionary effect of 
the cosmic epoch. Here we investigate the dependence of FR I and II dichotomy of 
radio galaxies on luminosity and redshift by using the SCR sample, where the FR I 
and II dichotomy was first seen, supplemented by data from two additional samples 
(MRC and B3-VLA), which go about a factor of 5 or more deeper in flux-density than 
the original SCR sample. This lets us compare sources with similar luminosities but at 
different redshifts as well as examine sources at similar redshifts but with different lu- 
minosities, thereby allowing us a successful separation of the otherwise two intricately 
entangled effects. 

Key words: galaxies: active - galaxies: nuclei - galaxies: evolution - galaxies: funda- 
mental parameters - radio continuum: galaxies 



1 INTRODUCTION 

One of the robust correlations in observational astronomy is 
between the morphology type of radio galaxies and their ra- 
dio luminosity. First pointed out by Fanaroff & Riley (1974) 
that there is a very sharp dependence of the morphology 
type of radio galaxies on their luminosity so that almost 
all radio galaxies below a luminosity Pits = 2 x 10^^ W 
Hz~^ sr~^ (for Hubble constant Hq = 50km s~^Mpc~^), 
are edge-darkened (called type I) in their brightness distri- 
bution, while all radio galaxies above this luminosity limit 
are more or less edge-brightened (called type II). This corre- 
lation has withstood the test of time (Miley 1980; Antonucci 
1993, 2012; Urry k Padovani 1995; Kembhavi k Narlikar 
1999). However, because of a very strong redshift-luminosity 
correlation (Malmquist bias) in a flux-limited sample, like in 



the SCR sample used by Fanaroff & Riley (1974), any effect 
related with redshift could appear as a luminosity-dependent 
effect and vice versa. One can then beg a question - could all 
differences seen in the FR I and II types of sources rather be 
a result of a cosmological evolutionary effect with the cosmic 
epoch, instead of, as almost universally believed, an effect 
of transition across a certain critical luminosity value? 

Following the archetypal paper by Fanaroff & Riley 
(1974), where they first pointed out the presence of two 
distinct morphology types of radio galaxies in the strong 
source SCR sample, ascribing the distinct morphology of 
each galaxy to its radio luminosity, it has ever since been 
thought to be only a luminosity-dependent effect (see e.g., 
Saripalli 2012 and the references therein). There have been 
no attempts to investigate the alternative possibility that 
it could as well be a redshift-dependent effect, demonstrat- 
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Figure 1. A scatter plot of both types of FR morphologies ('A' 
for edge-darkened FR I type and '+' for edge brightened FR II 
type) for the 3CR source sample of Fanaroff & Riley (1974) in 
the redshift-luminosity plane. To avoid undue compression of the 
scale, we have restricted the plot to a redshift limit of 0.15. There 
are no FR I type radio galaxies which lie beyond a redshift Of 
0.032 or have a luminosity Pits > 2 x 10^^ W Hz-^ sr-^ (for 
Hq = 50km s~^ Mpc~^) in the sample of Fanaroff &; Riley. 



ing the strong evolution of source morphology with cosmic 
epoch. For instance, suppose one were to correlate the mor- 
phology type of radio galaxies with redshift. Figure 1 shows 
a scatter plot with redshift and luminosity of both types of 
morphologies classified by Fanaroff & Riley (1974) in their 
sample. The demarcation is as good with redshift as it is 
with luminosity. It might be a moot point to guess what 
would have been the verdict, had Fanaroff & Riley in their 
seminal paper tried a correlation against only redshift, in- 
stead of luminosity. It is quite likely the effect then would 
have been interpreted as due to a very strong evolution with 
cosmic epoch, and the subsequent theoretical interpretations 
in that case perhaps very different. 

The two scenarios have very different physical interpre- 
tations, and either could be of equal importance. For in- 
stance in the conventional interpretation, with intrinsic lu- 
minosities being the root cause of their different morphology 
types, there is a huge amount of literature about the relation 
between the luminosity break and the different morphology 
types (Saripalli 2012 and the references therein). In fact all 
the models and discussions in the literature currently avail- 
able are almost exclusively only within that framework. On 
the other hand a definite correlation with redshift alone, 
would imply that it is the cosmic evolution of the properties 
of sources (perhaps because of the ambient density falling 
below a certain critical value due to the Hubble expansion) 
that might give rise to these two different type of morpholo- 
gies, with FR I type being the current favourite and the FR 
II type "a thing of past". All difference seen in FR I and 
FR II type sources, which are usually attributed to the dif- 
ferences in their luminosities, could as well be then related 
to the changing ambient densities with cosmic epoch. Also 
the very sharp division in luminosity could possibly be due 
to a critical ambient density value, which might divide the 
sources into two distinct morphology types. It might be the 
luminosity dependence or it might be the dependence on 



redshift that gives rise to these morphological differences, 
but this question could not be decisively settled based on 
any amount of arguments, sans actual observational data. 
At least this particular aspect has not yet been investigated 
in the literature. 

For one thing, to investigate this further one needs data 
from different samples to disentangle the luminosity and red- 
shift correlation in order to study the dependences of FR I 
and FR II break separately on them. That is to make these 
investigations one needs sufficient data comprising sources 
at different flux-density levels, so that one could examine 
sources with similar luminosities at different redshifts as well 
as compare them at similar redshifts for different luminosi- 
ties, thereby separating the two effects. We investigate this 
dependence of FR dichotomy of radio galaxies on luminosity 
and redshift by taking data from three such different sam- 
ples. We may add that there are reports of FR I types seen 
at redshifts larger than z > 0.5 (Saripalli et al. 2012), but 
a systematic investigation of this question is still needed us- 
ing samples which are complete in the sense that all FR I's 
above the sensitivity limit of the sample are included. 



2 THE SAMPLES AND THE DATA 

Our first sample is the SCR (Laing, Riley & Longair 1983), 
which is a complete strong source sample, with all necessary 
optical and radio information with good resolution maps so 
that one can in most cases unambiguously decide the FR 
I/II type of morphologies. 

The second sample we have chosen is essentially com- 
plete MRC sample (Kapahi et al. 1998) with ^os > 0.95 
Jy, which is about a factor ^ 5 deeper than the 3CR sample 
and has the required radio and optical information. Their 
total sample comprises 550 sources, with 111 of them be- 
ing quasars and the remaining are radio galaxies. Optical 
identifications for the latter are complete up to a red mag- 
nitude of 24 or a i^^ magnitude of ^ 19 and we do not 
think among the remaining there would be any missing FR 
I types, which in any case are not expected to be seen in 
this sample at redshifts larger than ^ 1. 

For our third sample we have chosen a homogeneous 
sub-sample, known as B3-VLA, from B3 sources with 5*408 > 
0.1 Jy, for which the relevant radio and optical information 
is given by Vigotti et al. (1989). Again the B3-VLA sample is 
expected to be complete at least as far as FR I type sources 
are concerned. Their VLA maps at 1415 MHz with 15 arcsec 
resolution are in general good enough for the morphologies 
to be classified into FR I/II types. 



3 RESULTS AND DISCUSSION 

To quantitatively distinguish between FR I and II, following 
Fanaroff & Riley (1974), we classify a radio galaxy as FR I 
if the separation between the points of peak intensity in the 
two lobes is smaller than half the largest size of the source. 
Similarly FR II is the one in which the separation between 
the points of peak intensity in the two lobes is greater than 
half the largest size of the source. This is equivalent to hav- 
ing the "hot spots" nearer to (FR I) or further away from 
(FR II) the central optical galaxy than the regions of diffuse 
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Figure 2. Histograms showing distributions of FR I/II morphol- 
ogy types radio galaxies with redshift for the (a) 3CR, (b) MRC 
and (c) B3-VLA samples. It should be noted that all radio galax- 
ies seen beyond the redshift limit of the plots (i.e. z > 0.25) are 
only of FR II type in all of the above three samples. 



Figure 3. Histograms showing distributions of FR I/II morphol- 
ogy types radio galaxies with radio luminosity (Pirs in units of 
W Hz-i sr-i) for the (a) 3CR, (b) MRC and (c) B3-VLA sam- 
ples. It should be noted that all radio galaxies seen beyond the 
luminosity limit of the plots (i.e. Pits > 10^^ W }iz~^ sr"-*^) are 
only of FR II type in all of the above three samples. 



radio emission. In all the three samples we have examined 
the radio maps and using the above criteria, we have classi- 
fied each source into either of the two types (FR I and II). 
There were a small number of sources that had a doubtful 
classification, we have dropped them from our analysis and 
that should not be too detrimental to our conclusions. 

Figure 2 shows histograms of the distributions of radio 
galaxies with FR I type morphology against redshift for all 
our three samples. We have restricted our plots to 2; = 0.25 
only, as no FR I type radio galaxy is seen beyond this red- 
shift in any of the three samples. We have determined me- 
dian value Zmed for the redshift distribution of FR I sources 
for each sample. The details of the method for determining 
the median values and for the estimation of their rms errors 
are described in Singal (1988). 

Figure 3 shows similar histograms of the distribution of 
radio galaxies with FR I type morphology against luminosity 
(at 178 MHz) for our three samples. Luminosity is calculated 
from their flux density 5*1 78 and the spectral index a {S (x 
y-"^) as 

Pl78 = ^178P'(l + ^)'+", (1) 

where V is the comoving cosmological distance calculated 
from the cosmological redshift z of the source. In general 



it is not possible to express T> in terms of ^ in a close- 
form analytical expression and one may have to evaluate 
it numerically. For example, in the flat universe models 
(Qm + = 1,^A 7^ 0), P is given by (see e.g., Weinberg 
2008), 



For a given Qa, V can be evaluated from Eq. (2) by a nu- 
merical integration. Here Qrn is the matter energy density 
(including that of the dark matter) and Qa is the vacuum 
energy (dark energy!) density, both defined in terms of the 
critical energy density = 3H§c^/{87tG), where G is the 
gravitational constant and c is the velocity of light in vac- 
uum. We have used Hq = 71km s"^ Mpc"\ Qm 0.27 
and Qa = 0.73 (Spergel et al. 2003). Since all the FR Fs 
we are interested in are at low redshifts {z < 0.25), the 
specific cosmological parameters do not make much differ- 
ence for the luminosity evaluations except for the scaling 
factor due to Hubble constant, which in our case makes lu- 
minosity estimates lower by a factor of (71/50)^ ~ 2 than 
the calculations in Fanaroff & Riley (1974). Thus the FR 
I/II luminosity break of Pits = 2 x 10^^ W Hz"^ sr"^ (for 
Ho = 50km s~^Mpc~^), arrived at by Fanaroff & Riley 
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Table 1. Percentile values of redshift and luminosity distributions for the FR I sources in the three samples. 

Sample 6'i78,med Number Z^ed ^Iq ^uq log(Pi78,med) log(Pl78,lq) log(Pl78,uq) 

Jy FRI 

SCR 19.4 23 0.03 lb 0.005 0.02 0.08 24.6 ib 0.2 24.2 25.5 

MRC 2.7 27 0.12 ±0.015 0.06 0.17 24.9 ±0.1 24.6 25.2 

B3-VLA 1.7 26 0.13 ±0.015 0.06 0.16 24.7 ±0.2 24.2 25.0 



(1974) in their seminal paper, corresponds to Pi 73 = 10^^ 
W Hz~^ sr~^ in our case. Again we have restricted our plots 
to Pi 78 = 10^^ W Hz"^ sr"^ only, as no FR I type radio 
galaxy is seen beyond this luminosity value in any of the 
three samples. 

From Figure 3, it is also clear that the FR I/II break is 
not as sharp as stated in Fanaroff & Riley (1974), as some 
overlap of both type of morphologies is seen in luminosity. 
However there is no denying that in all samples there are 
only a few, if any, FR I types with luminosities P178 > 10^^ 
W Hz-^ si-\ In Table 1 we have listed the median values 
both for the redshift and the luminosity in each of the three 
samples. Also listed are the median flux-density of FR I's 
in each sample and the number of FR I sources in the sam- 
ple. The rms error in Zmed in each case is determined from 
the frequency distribution (histogram) of ^-distribution. The 
rms error is given by y^/(2/p) in units of the class interval 
of z (Kendall 1945; Yule & KendaU 1950), where n is the 
total number of sources in the sample and fp is the ordinate 
value of the smoothened frequency distribution at the me- 
dian value in Figure 2. The median value of the distribution 
does seem to shift substantially with redshift for samples 
which differ in flux-density by a factor of 5 to 10. Pmed 
and the rms error in Pmed in each case is determined from 
the frequency distribution (histogram) of P-distribution in 
Figure 3, in the same way as for Zmed as described above. 
Within errors there is hardly any difference in three samples 
in the Pmed value, which thus seems to be independent of 
the flux-density level of the sample. 

While the median value of redshift {zmed) differs for the 
weaker samples as much as a factor of about four, at about 
a 5cr level, as compared to that in the stronger SCR sample, 
the difference in the median value of luminosity (Pmed) is 
only marginal, being only a factor of about 10°'^ ~ 2 and 
that too within only about la. From Zmed and log(Pmed) 
values, it is clear that the FR I type of morphology of radio 
galaxies is indeed due to their luminosity below a critical 
value as indeed envisaged first time by Fanaroff & Riley 
(1974) and that it is not directly related to the redshift and 
hence not due to a cosmic evolution effect. 

To ascertain it further we have examined the normal- 
ized cumulative distributions of FR I morphology type radio 
galaxies with redshift as well as luminosity. Figure 4 shows 
the normalized cumulative distributions of FR I morphol- 
ogy type radio galaxies with redshift for the three samples. 
From this figure we can compare the relative numbers of FR 
I types upto any given redshift for different samples, and 
thereby for different radio luminosities. If the morphology 
type depended only on redshift and not on luminosity, then 
all the three normalized cumulative distributions in Figure 4 
should be more or less coincident, which definitely is not the 
case. In Table 1 we have listed, along with the median value 
(^med), also the lower quartile (^^iq) and upper quartile (^uq) 
of the cumulative distribution of the redshifts of the sources 
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Figure 4. Normalized cumulative distributions of FR I morphol- 
ogy type radio galaxies with redshift for the (a) 3CR sample (con- 
tinuous curve), (b) MRC sample (broken curves) and (c) B3-VLA 
sample (dashed curve). Nl, N2 and N3 give the number of FR I 
radio galaxies in the SCR, MRC and B3-VLA samples, respec- 
tively. 



in the three samples. We see that the space distribution of 
FR I morphology type radio galaxies does change with the 
flux-density levels of the sample. There is a large difference 
in the redshift distribution of the SCR sample from those 
of the MCR or B3-VLA samples, the latter two (MRC and 
B3-VLA samples) however, have statistically almost indis- 
tinguishable redshift distribution of the FR I radio galaxies, 
not surprising as they are at about similar flux-density lev- 
els, while the 3CR sample is about 5 to 10 times stronger in 
flux-density than the other two samples. It does seem that 
the redshift distribution of the FR I radio galaxies depends 
heavily on the flux- density level of the sample. 

Figure 5 shows the normalized cumulative distribution 
of radio luminosity for FR I morphology type radio galaxies 
for our three samples. The three distributions overlap and 
there seems to be no gross difference in the luminosity dis- 
tribution of FR I sources in the three samples. In Table 1 we 
have also listed the lower quartile (Piq) and upper quartile 
(Puq) of the cumulative distribution of the radio luminosity 
of the source in the three samples, which again are of very 
similar values for all three samples. From these we find that 
half of FR I types lie in the narrow range of luminosities 
Pi78 ^ 1024-7±o.5 ^ jj^-i g^-i^ ^^^^ .g within a factor of 3 

around 5 x 10^^ W Hz~^ sr~^. Also all galaxies with radio 
luminosity above Pits ^ 10^^'^ W Hz~^ sr~^ are only of 
type II with an edge brightened morphology. We have thus 
reaffirmed that the FR I/II dichotomy is due to a change in 
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Figure 5. Normalized cumulative distributions of FR I morphol- 
ogy type radio galaxies with radio luminosity (Pits in units of W 
Hz"-*^ sr"-*^) for the (a) 3CR sample (continuous curve), (b) MRC 
sample (broken curves) and (c) B3-VLA sample (dashed curve). 
Nl, N2 and N3 give the number of FR I radio galaxies in the 
3CR, MRC and B3-VLA samples, respectively. 
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luminosity below and above a certain critical level, as first 
proposed by Fanaroff & Riley (1974) and followed in the 
literature ever since, and that it is not due to any cosmic 
epoch dependent evolution that gives rise to different FR I 
and II type morphologies. 



4 CONCLUSIONS 

We have investigated if there is a direct dependence of the 
FR I and II morphology types of the radio galaxies on red- 
shift. For this we compared their distributions in three differ- 
ent samples with different flux- density limits, which allowed 
us to separately study the effects of redshift and/or radio 
luminosity on the occurrences of the two morphology types. 
It was shown that the morphology type is not directly re- 
lated to redshift and thereby not a cosmic epoch dependent 
effect. The break between the two types of morphologies 
seems to depend only upon the radio luminosity. Half of the 
FR I type radio galaxies lie in the narrow range of luminosi- 
ties Pi 78 io24-7±o.5 ^ jj^-i g^-i^ ^.^j^ ^Q^g exceeding 
the value Pits ^ IQ^s-s ^ j^^-i ^^-i ^ above which all were 
found to be exclusively type II with an edge brightened mor- 
phology. 
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